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Abstract 

We explore the possibility of realising self-accelerated expansion of the Universe taking into account the vector 
components of a massive graviton. The effective action in the decoupling limit contains an infinite number 
of terms, once the vector degrees of freedom are included. These can be re-summed in physically interesting 
situations, which result in non-polynomial couplings between the scalar and vector modes. We show there are 
self-accelerating background solutions for this effective action, with the possibility of having a non-trivial profile 
for the vector fields. We then study fluctuations around these solutions and show that there is always a ghost, 
if a background vector field is present. When the background vector field is switched off, the ghost can be 
avoided, at the price of entering into a strong coupling regime, in which the vector fluctuations have vanishing 
kinetic terms. Finally we show that the inclusion of a bare cosmological constant does not change the previous 
conclusions and it does not lead to a ghost mode in the absence of a background vector field. 



1 Introduction 

Understanding the nature of dark energy represents one of the most interesting open questions in cosmology. Mod- 
ifications of Einstein's theory of gravity at large scales could explain the present day acceleration, without invoking 
the aid of a cosmological constant or a exotic matter content. Therefore, self-interactions in the gravitational sector 
may be sufficient to realise cosmological acceleration with no need of an additional energy momentum tensor; this 
phenomenon is dubbed self-acceleration. 

Massive gravity is an example of these modified gravity models, in which the self- acceleration is realized thanks 
to the specific dynamics of gravitational degrees of freedom. Already in the Fierz-Pauli version of massive gravity 
[1] , explicit solutions were found describing a de Sitter spacetime [2] , with the de Sitter radius inversely proportional 
to the graviton's mass. Recently, in the context of a novel non-linear extension of Fierz-Pauli [3, 4], self-accelerating 
solutions have been discovered [5, 6, 7, 8, 9, 10] [See also [11, 12, 13, 14, 15] for a non-exhaustive list of work in this 
theory and related models]. In a healthy theory of massive gravity, the graviton contains five degrees of freedom: 
two tensors, two vectors, and one scalar. The dynamics of the scalar mode has been studied in detail in [8], 
showing that a non-trivial configuration for this field leads to self-acceleration. The analysis was carried out in 
the so-called decoupling limit, a regime of linearised dynamics in the tensor modes, and in which the scalar-scalar, 
as well as scalar-tensor interactions, are described by a finite number of terms, which can be transformed into 
the so-called Galileon combinations [16]. Scalar fluctuations around these self-accelerating configurations are free 
of ghosts, provided that parameters characterizing the theory are chosen appropriately [8]. On the other hand, 
around this self-accelerating configuration, the vector degrees of freedom enters into a regime of strong coupling, 
as the coefficient in front of the kinetic term vanishes [8]. 



However, in all the literature of self-accelerating solutions in massive gravity, little attention has been put into 
the vector mode. One may ask if the theory admits other classes of self-accelerating configurations, which include 
vectors, and if they exist what is their perturbative stability. In this paper we answer these questions, presenting 
new solutions in which the vector degrees of freedom play a crucial role, supporting, together with the scalar mode, 
the self-acceleration. We focus on the theory in the decoupling limit, and develop tools to obtain an effective action 
describing the vector modes in this limit. This effective action turns out to have an infinite number of terms for 
the vector-scalar interactions, in contrast to the finite number for the non- vectorial sector. This infinite series leads 
to non-polynomial interactions between the scalar and the vector. By concentrating our attention in spherically 
symmetric ansatze, we are able to resum the series, and find general self-accelerating solutions in the effective theory 
which describe de Sitter space with non-trivial profiles for the scalar and vectors. These general solutions, as far 
as we know, represent the first example of self-accelerating configurations with a vector field. We then proceed to 
study the dynamics of fiuctuations around these general solutions. We find that a non-trivial profile of the vector 
field is able to alleviate the strong coupling behaviour of vector fluctuations, by providing positive definite kinetic 
terms to them. However we find that one of the scalar and vector perturbations always becomes a ghost, so that 
the only ghost-free self- accelerating solutions are those without the background vector field. Furthermore, a bare 
cosmological constant can be included in these self-accelerating solutions and the theory of fluctuations about them 
remains ghost-free if the background vector mode is not present. However, in the case of a non-trivial profile for 
the background vector field, the ghost mode cannot be avoided. 

The paper is organized as follows. In Section 2, we review the construction of the non-linear Lagrangian for 
massive gravity. In Section 3, we focus on the simplest version of this theory, presenting new self-accelerating 
solutions with a vector field included. In Section 4, we determine the Lagrangian for massive gravity in the 
decoupling limit including vector degrees of freedom. In Section 5, we study the dynamics of linear fluctuations 
around the self-accelerating configurations, showing that a ghost is always present in the spectrum if there is 
background vector. In Section 6, we extend our analysis to the full theory including two additional parameters 
and determine the most general self-accelerating solutions in the decoupling limit with vector fields. The analysis 
of linear fluctuations around these configurations show that there is no parameter space where we can remove the 
ghost if a background vector field is switched on. We present our conclusions in Section 7. Three Appendixes 
complete the paper with technical details. 



2 Massive Gravity Lagrangian 

We consider the following Lagrangian for massive gravity, a non-linear extension of Fierz-Pauli theory proposed in 
[4] 

£=^V^{R-2A-U). (2.1) 

The potential depends on a dimension-full parameter m, which sets the graviton mass scale, and on two dimen- 
sionless parameters as and a^. It has the following form: 



U = -m^ [U2 + azUs + ai U4] , (2.2) 



with 



U2 = (tr^)2-tr(/C2), 

Us = (tr/C)3-3(tr/C)(tr/C2) + 2tr/C^ 

U4 = (tr/C)''-6(tr/C)^(tr/C2)-h8(tr/C)(tr/C3) + 3(tr/C2)2 -etr/C''. 
The tensor K.^ is defined as [4] 

IC; EE S;-[V9-'[9-h])^^ , (2.3) 

/ ^ I ^ 

where the square root of a tensor is defined as vA^^ vA^^ = A^^'^, for any tensor A^^". The metric h^^ is 
the displacement from a fiducial flat metric r/^^, namely /i^^ = t\^i, — 5^^. The action (2.1) breaks explicitly 
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reparametrisation invariance, but it can be restored by means of the Stiickelberg trick [11]. The new definition for 
the tensor JC^ is given by 



IC,"" ^ 6;-{Vg-^g-H])^ , (2.4) 



where now iJ^^ corresponds to the covariantisation of the metric perturbations, defined as 

H^i, = hf,„ + i]p^dfj,n'^ + r]a^id^i^°' - 'rjafidf.n^'di.TT'^ . (2.5) 

Therefore, a change of coordinates a;^ + should be accompanied by the following transformation of the 

Stiickelberg field tt'^, 

vr'^^Tr'^ + e'', (2.6) 

in order to recover full diffeomorphism invariance. The field tt^ can be decomposed into a divergence-less vector 
A^, and a scalar tt in the following way 

tt'^ =,yA'-(A,+5,7r). (2.7) 
Consequently, in terms of the vector and scalar, the tensor iJ^jy can be expressed as 

H^,, - h^,. + 2n^, - + d,A. + d,A^ - d^A'^d.A^r^^p - n^„9,A" - 9^^"n,„, (2.8) 

where 11^^ = n^all" and H^^ — df^di^TT. The indices of tt^ are raised/lowered by the flat fiducial metric ry^ii/. 

The Stiickelberg scalar tt and vector A'^ play the roles of scalar and vector components of the massive graviton: 
the Lagrangian obtained by applying the Stiickelberg trick contains the same number of degrees of freedom as the 
original Lagrangian (2.1). By plugging the expression for the tensor /C^'', written in terms of h, A and tt, inside 
(2.1), one obtains a Lagrangian for the massive gravity theory expressed in terms of tensor, scalar, and vector 
degrees of freedom. In order to canonically normalize the degrees of freedom^, we need to rescale the fields as (see 
for example [12]) 

hf,^ Mpi hf,^ , mMpi Af, , tt -> m^Mpi tt, , A ^ Mpi A. (2.9) 

The transformations associated with the diffeomorphism invariance Eq. (2.6) can also be rewritten in terms of these 
fields as 

Mpi 
IVlpi 

Jtt = -mA, (2.10) 

for an arbitrary divergence-less vector and a scalar A. This Lagrangian contains various non-linear interactions 
characterized by derivative couplings, which in turn play a crucial role for allowing an implementation of the 
Vainshtein mechanism [5, 6, 17, 13], and for generating self-accelerating configurations leading to the de Sitter 
expansion [3, 5, 6, 8, 9, 10]. 

The aim of this paper is to discuss new self-accelerating solutions for massive gravity, exploiting the features of 
the vector sector of the theory. In order to do so, it is convenient to focus in the decoupling limit [11] 

TO ^ , Mpi oo , A3 EE m^Mpi = fixed. (2.11) 



^For the tensor h^,^ and vector A^, the kinetic terms set the canonical normaUzation, whereas for the scalar field tt the kinetic 
terms are total derivatives. Given that one would like to keep some terms involving the scalar field in the decoupling limit (2.11), 
we use normalization (2.9). This preserves non- vanishing scalar-tensor couplings in that limit, hence providing kinetic terms for the 
scalar after a diagonalisation. Similarly, the cosmological constant A should scale as hfj,,, in order to have consistent solutions in the 
decoupling theory. 
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This limit is particularly appropriate for our objectives, because the tensor piece in the Lagrangian simplifies 
considerably. On one hand, the tensor self-interactions are described by the quadratic expansion of Einstein- 
Hilbert action, while interactions with other degrees of freedom are linear in the tensor mode. On the other hand, 
the Lagrangian describing vector and scalar degrees of freedom maintains its non-linear structure, admitting de 
Sitter solutions with no need of an additional energy momentum tensor [3, 5, 6]. In these solutions, the de Sitter 
radius is proportional to the inverse of the strong coupling scale A3 . 

The Lagrangian in the decoupling limit of Eq. (2.11) has a reduced symmetry with respect to (2.10), namely 



hence the symmetry is reduced to the linearised diffeomorphism invariance times an independent U{1) symmetry 
for the divergence-less vector. The Lagrangian in the decoupling limit, when expanded from a flat fiducial metric, 
has been shown to be free from the so-called Boulware-Deser (BD) ghost [18, 3], in the sense that it contains only 
five degrees of freedom; the sixth mode, potentially leading to a ghost, is absent ^. However, around the self- 
accelerating de Sitter configurations, one of the remaining physical five degrees of freedom may become a ghost. 
As we will learn in this paper, there are conditions on the parameters of the theory in order to have a set-up that 
is ghost free around the self-accelerating configurations. 

There are strong indications that the theory is ghost free (in the sense that the Boulware-Deser sixth mode is 
absent) also away from the decoupling limit [14], although there is still debate on this issue [19, 20, 21, 22]. The 
analysis of this subject is not the main aim of this work, so we will not discuss it here. 

In the next section, we present a method for determining new self-accelerating configurations, in the decoupling 
limit of this massive gravity theory, but including the dynamics of vector degrees of freedom. We determine them 
in two ways; firstly, we show that they can be obtained by taking the decoupling limit of known exact solutions 
in the full theory. Secondly, by imposing the appropriate symmetries, we derive an effective action for the relevant 
degrees of freedom (including vectors) whose general solution provides the same self-accelerating configurations. 

3 Self-accelerating solutions with vectors {a^ = = case) 

In order to study new self-accelerating solutions in the presence of vector degrees of freedom, we proceed as follows: 
we start with the most general static spherically symmetric solutions for the original Lagrangian (2.1). As this 
Lagrangian stands, it describes only the dynamics of the metric h^j^^, and as we previously showed [5, 6], there is a 
branch of solutions exhibiting the static patch of the de Sitter spacetime. 

By reintroducing the diffeomorphism invariance via the Stiickelberg trick (2.4), we apply a gauge transformation 
that allows to recast the fields in a form that is suitable for taking the decoupling limit (2.11). The details on how 
this procedure is performed are in Appendix B. After this limit is taken, we get a self-accelerating configuration 
where the de Sitter invariance is fully manifest. Moreover, we find that these solutions include a non-trivial profile 
for the vector field. As far as we are aware, this is the first example of self-accelerating configurations which includes 
vector modes. 

We then proceed to show that the same class of self- accelerating configurations can be alternatively obtained 
as solutions for the Lagrangian in the decoupling limit, imposing an appropriate Ansatz for the fields involved that 
respects the symmetries of the de Sitter spacetime. 

3.1 Self-acceleration from spherically symmetric solutions 

The problem of classifying the most general spherically symmetric static solutions for the original, non-covariant, 
Lagrangian (2.1), describing the dynamics of tensor degrees of freedom, has been already discussed in the literature 

^When one considers the fiducial flat metric as the background for studying linear fluctuations, there is no way of exciting a possible 
sixth degree of freedom using the vector-scalar coupling, since the action in the decoupling limit is quadratic in the vector field. However, 
as we show in this paper, vector fluctuations play an important role around different backgrounds. 



6h 



0, 



Sir 



(2.12) 
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[5, 6, 23]. If de Sitter solutions exist, their static patches (if any) have to be contained within our classification. 
Here we briefly review our findings. We start with the most general Ansatz for static spherically symmetric 
configurations, 

ds'^ = -C(r) dt^ + A{r) dr^ + 2D{r) dtdr + B{r)dfl^, (3.1) 
where dfJ^ — dO"^ + sin^ OdcjP' . The non-dynamical flat metric is written in terms of spherical coordinates as 

ds)iat = -dt^ + dr'^ + r^dn^ ■ (3-2) 

Notice that in general relativity, one can set B{r) = by a coordinate transformation, but this is not possible 
here, since we do not have diffeomorphism invariance. In order to simplify our analysis, it is convenient to define 
the combination A(r) — A{r)C{r) + D^{r). We plug the previous metric into the Einstein equations 



Gp,iy — T^iin (3-3) 



where the energy momentum tensor associated with the potential U of Eq. (2.2) is defined as T^^^ — ^23— ^^^~s'^ 



The Einstein tensor G^^ satisfies the identity D{r) Gu + C{r) Gtr = 0, which implies the algebraic constraint 

= i?(r)r^ + c(r)r^ 

D{r) (2r - 3/BW) 
= m^^—^ ^ .1/. - (3-4) 

The previous condition can be satisfied in two ways, which lead to two different branches of solutions [5, 23] (see 
also [24, 25]). We can either set D{r) = 0, and focus on diagonal metrics, or alternatively set B{r) = 4r^/9. The 
fact that there are two branches of solutions indicates that, unlike in general relativity where the Birkhoff theorem 
holds, there is no uniqueness theorem for spherically symmetric configurations in this theory. The diagonal branch 
gives solutions which are asymptotically flat [5, 6]'^. Since we are interested in de Sitter configurations whose 
curvature invariants do not decay at infinity, we will not discuss further the diagonal branch of solutions. On the 
contrary, the second branch provides asymptotically de Sitter configurations. To show this explicitly, we impose 
B{r) = 4r^/9. Due to identity C{r)Tl^^ + A{r)T^f = 0, one gets a further condition 

A(r) = A{r)C{r) + D^{r) = Aq = const. (3.5) 

The remaining Einstein equations provide the following unique solution (see Ref. [6] for detailed derivations) 

A{r) = ^(p(r)+a + l), B{r) = ^r^ (3.6) 
Cir) = ^(l-p(r)), i?(r)^^AoVp(r)(p(r) + a), 

where 

with arbitrary constants c and Aq. Notice that this configuration depends on two integration constants. A sufficient 
condition to ensure that D{r) is real, is to choose c > and < ^/Ao < 4/9. The integration constant c corresponds 
to the Schwarzschild mass, but since we are interested in pure de Sitter metrics, we set c = in what follows. The 
general case (as, a4 ^ 0) can be analysed in the same way, but its discussion is postponed to Section 6. These 
solutions correspond to a maximally symmetric de Sitter space, since the Ricci tensor is proportional to the metric. 
If we perform coordinate transformations and go to the flat Friedman-Robertson- Walker slicing of the de Sitter 
spacetime, there appears a coordinate singularity at the de Sitter horizon. The physical nature of these singularities 



^The decoupling limit Lagrangian of this branch has only one asymptotic solution, which indeed decays to the fiat space at a large 
distance from a source for the case 03 = 04 = [5]. However, in the most general case, where 03 and 04 do not vanish, there are more 
asymptotic solutions than the flat space one [6]. 
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is an interesting open issue, however since we focus in the decouphng limit, this singularity does not show up in 
our analysis. In fact, it was shown that there are no closed or flat FRW solutions in this theory [9]. However, a 
self-accelerating open FRW solution was recently found [10], and in the decoupling limit, this last solution is indeed 
included in the solutions that we obtain from exact spherically symmetric solutions. 



As explained in the previous sections, we can restore the diffeomorphism invariance via the Stiickelberg trick, 
and apply a gauge transformation that recasts the metric of Eq. (3.1) (with c = 0) in a form that makes the de Sitter 
nature of this space-time manifest. We developed this procedure in [5, 6], where we wrote the de Sitter spacetime in 
a explicitly time-dependent frame, while in what follows we will express it in conformally flat coordinates'* . Details 
on how to obtain the decoupling limit of solution (3.6) can be found in Appendix B, but here we only describe the 
final result: 



ds^ = 
A, = 

TT — 



1 - 



'PI 



{-dt^ + dr"^ + r'^dVt^) + 0{m^) , 



Aar^ /16 



24 V A, 



A, 



3 ..2 ^2. 



81 + 0(771^), 0,0,0 , 



3 A, 



1 t^ + Oim'). 



(3.8) 



where we used Mpi = h.^/m? and only wrote the term at leading order in m, since we intend to focus in the 
decoupling limit (2.11). Indeed in this particular limit, we take m — >■ 0, so only the terms we explicitly wrote 
survive. The metric then assumes a manifestly de Sitter form in a conformally flat frame. This configuration 
provides a generalization of the solution analysed in [8] where only the scalar field is switched on. Actually, the 
choice Aq = 16/81 switches off the vector and the solution reduces precisely to that of [8]. However, in general 
also the vector degree of freedom is present. Moreover, it is important to stress that even though we derived these 
solutions (3.9) from the specific non-linear solution (3.6), the decoupling limit solutions are more general in the 
sense that many different non- linear solutions can be reduced to the same decoupling solutions (3.8). For example 
and as mentioned before, the open FRW solution obtained in Ref. [10] reduces to (3.8) with Aq = 16/81 in the 
decoupling limit. We will study various full non-linear solutions describing the self-accelerating universe in a future 
publication [28]. 



3.2 Solving the field equations in the decoupling limit 



It is also possible to reproduce the same configurations by solving the equations of motion for the Lagrangian in 
the decoupling limit by imposing the appropriate symmetries. In order to maintain spherical symmetry, we adopt 
the following Ansatz for the canonically normalized fields. 



A,, 



1 



Mpi 
(A(r), 0,0,0), 

TTt{t) + TTr{r)- 



(3.9) 



This Ansatz is able to generate the de Sitter metric expressed in conformally fiat coordinates. Although it is not 
the most general Ansatz compatible with the required symmetries, it is sufficient for reproducing the solution (3.8). 

After inserting the Ansatz (3.9) in the action (2.1), we find the following Lagrangian in the decoupling limit for 
the fields /, At and tt (with the fiducial flat metric expressed as (3.2)), 



-4'^(/) 



-F2in)f{t,r) + F3in){drAt)' 



(3.10) 



*Note that these coordinate transformations can unnecessarily introduce singularities in tt^ 
singularities do not show up. 



Again in the decoupling limit these 
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where 



W = 2(3+^)^+(3+i^)."-(3 + 2."+4)., 

The dot and the prime correspond to derivatives along t and r respectively. Moreover, we set A3 = 1 for convenience. 
The r^Fi/m? term is a total derivative, thus it does not contribute to the equations of motion. The second term, 
Cj^ , is the Einstein-Hilbert action truncated to quadratic order in /. The direct coupling between / and tt (i.e. 

F2(t^) f) corresponds to the contribution found in [3], and denoted there by h^'^ {x\^} + X^J), where X^""* has n 
powers of df^d^TT. This coupling leads to the second order differential equation for the field tt. 

Regarding the vector sector, we obtain a non-polynomial coupling between the vector and scalar degrees of 
freedom controlled by the function i^3. Derivatives of the scalar field appear in the denominator of this expression 
even in the decoupling limit. We will see that this is a result of a resummation of an infinite number of terms. We 
will discuss how to obtain these couplings in general in Section 4. The equation of motion for the scalar mode, tt, 
is given by 

The first term in the previous expression could lead to higher order derivatives in time if SC/Stt depends on tt. The 
F2 coupling in the Lagrangian will not generate such higher derivatives, in agreement with the findings of [3]. In 
contrast, the vector-scalar coupling might give rise to such higher derivative terms, which might be thought as 
the propagation of the sixth-degree of freedom. However, as we will see later, in the equations of motion describing 
fluctuations around our self-accelerating solutions all these higher derivative terms cancel exactly, ensuring that 
the Boulware-Deser mode is absent. 

The equations of motion associated with Lagrangian (3.10) have the following asymptotically non-decaying 
solution 

A3 



A 



A3 Qo 2 



2 

3 
4 




where Qo > is a free integration constant. It is straightforward to check that, choosing Qo — y ^ " the 

previous solution indeed reduces to (3.8) when m 0. This proves that the configuration (3.8) is a solution of the 
field equations in the decoupling limit. 



4 The Lagrangian with vector degrees of freedom 

We are now interested in deriving a more general expression for the Lagrangian in the decoupling limit including 
vector degrees of freedom. As we will learn, this Lagrangian contains an infinite number of terms, however, this 
infinite series of terms can be resummed for the self-accelerating solutions. 

4.1 The potential in the decoupHng Hmit 

In order to derive a Lagrangian for the decoupling limit of the action (2.1), we write Mpi in terms of m using 
(2.11), and then take the m — )• limit. To avoid to be overburdened with indexes, it is convenient to write tensors 
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in a matrix form; a given tensor S^j, may be written as B. Using this matrix notation, the tensor Hfj^^ = H reads 



H = g - {1 -n?]-^ - dAr]-^) {1-117]-^ - dAr]-^) , (4.1) 

where g — g^^, rj = rj^^, r]~^ — rj^'^, H = 11^^, dA = S^A^, and 1 is the identity matrix. 

The quantity we are interested in is IC^^, = g^^K. = 1 — -^/l — g^^ H; so it is more convenient to write 1 — g^^H 



as 



1-g-^H = g-'^ [I-Ut]-^ - dA'q-'^) {rj-U- dA'^ 



(4.2) 



= (1 - rn^h-q'^) |(1 - -qUf- - rmj^'^dA (1 - rjU) - m {1 - f]n) 'q-'^dA'^ + m'^ r]-^dAr]-'^dA'^^ , 

where in the last step we have written everything in terms of the canonically normahzed fields (2.9), and expand 
the expression up to second order in m. There are also factors of the strong coupling scale A3 defined in Eq. (2.11) 
that we set to one. However, it is straightforward to put it back by means of dimensional analysis. 
We define 



so that we can succinctly write 



p = (i-7?-in), 

Li = T]-^ dA P + P rj-^ dA^ , 
L2 = t]-^ dAi^-^dA^ 



(4.3) 
(4.4) 
(4.5) 

(4.6) 



We can neglect tensor fluctuations since, in the decoupling limit, they only couple to the scalar but not to the 
vector mode. The scalar-tensor couplings have been already completely classified in [3], so we can use these results 
to complete the full Lagrangian. Therefore, from now on we set h = Q. 

The potential for the scalar and vector fluctuations we are interested in is given by (2.2). Without loss of 
generality, we focus on the — a^^ = Q case first, and then generalise the result to non- vanishing 03 or in 
Appendix C. For = = 0, the potential (2.2) simply reads ((/C) = tr/C) 



U 



where, again, we set A3 = 1. Calling (l — 5 ^ i?) = M so that g ^/C = 1 — \/M, we obtain 



6(VM) + (M) - {VMy - 12 



(4.7) 



(4.8) 



Thus all the problem is to calculate trvM, to second order in to. This will be enough since the potential (4.7) 
has an overall factor of m^^ and contributions appearing with powers of to bigger than two will go to zero in the 
decoupling limit. Therefore, up to second order in m, we can write 



M = P' [1 - mP"2 Li + TO^P-' L2] 

^2 



P' 



m 



P2 [1. 



2 

mQi 



P-'L, 



-P-'Lo 



with 



Qi 



p-^Li 



(4.9) 

(4.10) 
(4.11) 
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Thus, the square root of M can be written, up to 0{m'^), as 

Vm = P [l — m Qi + rri^ (52] + m D + rri^E, (4-12) 

where the matrices D and E are needed since P and 1 — m Qi + Q2 do not commute; D and E obey the 
equations determined by requiring that the square of equation (4.12) reduces to (4.9), to second order in m. Using 
(4.12), one arrives to the foUowing expression for the potential U, 

U = Ati [PQ2 + E]+tTL2+0{m) (4.13) 
= -itr [P-^ rj-^dArj-^dA + p-^ rf^dA^ ij-^dA^ - p-^ r^-^dA-q-^dA^ + r]-^dA^ p-^ r^-^dAP] 

+tr [r^-^dAr]-'^dA^] + Air E + 0{m), (4.14) 

where we have neglected total derivatives and used that tr = 0, which we will prove later. Notice that all the 
quantities proportional to the inverse powers of m disappear, and when taking the decoupling limit the 0{m) term 
vanishes. However, we still need to determine tr E. As mentioned before, by taking the square on both sides of 
Eq. (4.12), we get the following condition to second order in m, 

= -mP[{Qi~mQ2),P]{l-mQi) + m{D,P{l-mQi)} 
+m^{P, E} + 771^0^. 

To first order in m, we find a matrix equation for D given by 

{D,P} = P[QuP], 

that can be recast as 

D + P-^DP = [Qi,P]. (4.17) 

The previous equation implies that tr D = 0, as promised. In terms of the original vector and scalar fields, equation 
(4.17) can be rewritten as 



(4.15) 
(4.16) 



D + PDP-^ = ^{p-^ir'^dAP-r]-^dA) - ^{PT]-^dA^ P-^ -Tj-^dA^ 



(4.18) 



In certain cases, in which P and A have particularly simple forms, a solution to the previous matrix equation can 
be guessed by inspection of the structure of the matrices. On the other hand, it is convenient to have a systematic 
method to deal with solutions of the previous equation. It is easy to check that a formal, general solution can be 
written in terms of the following combinations of infinite series 



2D = ^a„F-"77-idAP" + ^6„P"7;-idAF-" 

n=0 n=0 
CO 00 



(4.19) 



T!=0 ri=0 

where the coefficients in the series satisfy the following relations 

a„+i = for n > 2 , 



bn + bn-i — for n > 2 , 
ao + ^0 + ^1 = 
ao + Bq + ai = -1 

fli + a2 = 1 



and 



c„ + c„+i =0 for n > 2 , 
dn + dn-i = for n > 2 , 
co + ci+do = 1 
co + do + di = -1 

Cl + C2 =0 



Notice that the previous relations can be satisfied in different ways. For example, setting the coefficients 6„ = 
dn, one obtains 



2D 



-P-^r]-'^dAP + 2 (-1)" p-^rj-'^dAP'' 

n=2 

00 

+ rj-^dA^ + 2 (-1)" ^" v'^dA^ P"" . 



(4.20) 
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We will later discuss a method that allows to resum these series, but for the moment, let us continue determining 
the quantity trE, that is needed to calculate the potential in Eq. (4.14). From Eq. (4.15), to second order in m, 
one finds 

= P [Q2, P]+P [Qi, P]Qi-{D,PQ,} + {P, E} + D^, (4.21) 
which after some manipulations and using (4.17), (4.21) can be put into the form 

= [g2, P] + [D, Qi] + P-^ {P, E} + D^. (4.22) 

Taking the trace, one finds 

tiE = -itr {P-^ D^) . (4.23) 

In order to get tr£', the only unknown quantity'' in (4.14), one needs to determine tr (^P~^ D^) . Plugging the 
solution for D of Eq. (4.20), the previous quantity involves the resummation of an infinite series. Therefore, we 
learned a very interesting fact: even in the decoupling limit, the Lagrangian describing vector degrees of freedom 
contains an infinite number of terms. This has to be compared with the case of scalars only, in which we only have 
a finite number of terms corresponding to the Galileon combinations [3] . 



4.2 A method for resuming the series 

We now discuss a method, valid in many physically interesting cases, for resuming the infinite series we encountered 
in the previous subsection. This method aims to rewrite those series in terms of geometric series. These can be 
easily resumed, leading to the non-polynomial vector-scalar couplings that we met in Section 3.2. 

We make the hypothesis that the matrix P is diagonalisable. We write it as 



P = UTU' 



(4.24) 



where J7 is a unitary matrix, while T is a diagonal matrix that has the eigenvalues Vi of P in the diagonal. 
Eigenvalues and eigenvectors of P are particularly simple to calculate in the spherically symmetric case, since the 
matrix P splits into 2x2 matrices with particularly simple structure. For more general cases, extra work is needed. 
As seen in (4.19), we would like to calculate the quantities 

p-"r]-^dAP", P"Tj-^dAP-'\ P~'^r]-^dA'^ P", P"" t]'^ dA^ P"\ (4.25) 

Considering for example the first of the previous quantities, we can write 

p-'^rj-^dAP" ^ [/T-"BT"C/-\ (4.26) 

with 

B = U-^r,-^dAU = {&pJm' 

where {bpq}J^.J indicates a matrix with a component bpq at the p—th row and the q-th column. Then, we use the 
ity 

k' 



following identity 



'n-l 

E 

.k=0 



- 1 



(4.27) 



The case n = 1 is easy to check, and for general n, it can be proved by induction. The geometric series in the 
previous formula can be then resumed, thus one gets 



T"" ST" = B 



1 



M 



(4.28) 
(4.29) 



M 



^As we will see in section 6, in case where 03 or 04 do not vanish, we need the complete expression for E, and not only its trace, 
since the potential depends on M'^/'^ as well as \/M. 
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This is what we need to resum exphcitly the series appearing in the expression (4.19), which reduces to a combination 
of geometric series. 

To conclude, once the eigenvalues and eigenvectors of P are known, our method allows to rewrite the infinite 
series contained in tr E in terms of geometric series that can be straightforwardly resumed. Instead of providing 
general but lengthy formulae, we will discuss concrete and interesting examples in what follows. 



4.3 An example 



Let us apply the previous method to a representative example. We adopt spherical symmetry, and write a particular 
Ansatz which describes the configuration as in (2.10). Thus, we consider the following setup 



= (At(r), 0,0,0), 

TT = TTt{t) + TTrir). 



Then 



P = 



1 + 7f 








' \ 





1 -tt" 














1 - Tr'/r 














1- Tr'/r / 



/ 



V 



-A'tir) \ 





/ 



(4.30) 



(4.31) 



Using the formulae introduced in the previous section, it is straightforward to check that 



p-^rj-^dAP'' 
P" ri-^dA'^ P-" 



1 



Tr 



1 + TT 

1 - Tr" 



1 



■ Tr 



r]-^dA, 
Tj-'dA^. 



(4.32) 
(4.33) 



Applying these results to the solution for D written in Eq. (4.20), we can easily see that those series become the 
geometric series. We can indeed write 



2D 



-P-^ rj-^dA P + ri-^dA^ + 2 irj-^dA + rj-^dA^ 



n=2 ^ ' 



(4.34) 



The series is convergent if |(1 — Tr")/(1 + Tf)| < 1, a condition that, a posteriori, can be checked to be satisfied 
for our selfaccelerating configuration (2.10). Using the usual formula for resuming the geometric series, we get the 
explicit solution for D as 



(Tr" - 1) (ii- + Tr") 



2(l+ii-)(2 + ii--Tr") 



2(2 + ^- Tr") 



(4.35) 



Substituting the solution for D in the expression for the trace of E, Eq. (4.23), one recovers the scalar- vector 
coupling in the Lagrangian (3.10). The appearance of scalar derivatives at the denominator of the expression is 
now easy to understand; it is due to the resummation of an infinite number of terms in the geometric series. 

Before moving on to the stability of the self-accelerating solution with a vector charge, we consider at this stage 
a more general Ansatz preserving the spherical symmetry, which will be necessary for the perturbation analysis. 
Consider 



Tr(i, r) 



and 



{At{t,r),Ar{t,r), 0,0) 



as Ansatz for the scalar and vector that preserves the spherical symmetry. Then the resulting interaction term 
between the scalar and the vector, calculated using the previous resummation techniques, is given by the following 
contribution to the action 



Ssv = 



r'^dtdr< 



■2{l + n- Tr" - 2^)ArA', + (1 + 2^) (i^ + 
2(2 - Tr" + 



- A'^At + -Ar{K- At) + ^Al}, (4.36) 



1 
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which after integration by parts and using the divergence-free condition for the vector ((r^A^)' = r'^At), takes the 
simpler form 

1 f r^dtdr / 27r'\ / • \2 



5 Stability of the self-accelerating configuration (0^3 = a4 = case) 

In this section, we analyse the perturbative stability of the self-accelerating configm^ation (3.8), in the simple case of 
as = a4 = 0; we will discuss the most general case in the following section. For simplicity, we consider perturbations 
that preserve the spherical symmetry. We discuss their dynamics, and investigate the existence of ghost degree 
of freedoms around our configurations. We separate the discussions of this section in different subsections, which 
focus on different contributions to the total Lagrangian. In the last subsection we collect the results and study the 
sign of the kinetic terms for the relevant modes. We show that a ghost is always present around the self-accelerating 
configuration with a background vector field. 



5.1 Scalar- Vector Lagrangian 

We can decompose the divergence-free vector in terms of scalar and vector components with respect to the 
2-sphere. Vector components are free and do not couple to tt, so we can neglect their dynamics while studying the 
stability of the background. Thus in the following we only consider the scalar perturbations with respect to the 
2-sphere. The residual U(l) gauge freedom of Eq. (2.12) allows us to choose a gauge for so that it has only time 
and radial components. As a result, the most general Ansatz for scalar and vector field that preserves spherical 
symmetry for our system is 

TT = 7r(i, r), A^^ (At(t, r), A,(t, r), 0, 0). (5.1) 

The action controlling interactions among the scalar tt and vector A^ in the decoupling limit, is given by (4.37). We 
now consider the spherically symmetric, time-dependent perturbations around our self-accelerating configuration 
(3.12), namely 

TT^^o + M^,^), A^ = (Ato(r) + Mt(i,r),M,(i,r),0,0). (5.2) 



where the background is given by 




TTo = -\[r^-e)+^^\\ll+^-^-l]t^, (5.3) 
A„ = f (5.4) 



The action (4.37), expanded at second order in perturbations, can be split as 



Ssv = 4y + S'^sv (5-5) 



where 



4y = / dtdrr' 



9 + 16Q§)(57r' 



(5.6) 



is linear in the perturbations, and 



c(2) ^ qS I qint (r rj\ 

^SV — '^SV + "^SV ^ 10- 'J 



oS 



Qint 



dtdrr 



dtdrr^ 



(3 + v/9 + 16Q2): 

:(57r'(M; - 5Ar) 



3 -f V9 + 16Qo 



(5., 



(5.9) 
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is quadratic in the fluctuations. Using the divergent free condition of the vector 

r 

we can show that Sgy can be re-expressed as 



?(2) 



/ 



dtdrr 



4Qo 



(3 + v/9+16g2)2 



-^QoiST:' + (Jtt^) + (3 + -y/g + 16Ql){67r'SA't - bixbAt) 



(5.10) 



(5.11) 



As noticed before [11], if no background vector field is present there are no dynamical vector perturbations at the 
quadratic level. However, any small non- vanishing Qqi produces dynamics in the vector degrees of freedom. Notice 
also that a non-vanishing Qq spontaneously breaks Lorentz invariance of the perturbed action in the scalar sector. 



5.2 Scalar-tensor Lagrangian 

In the decoupling limit, the scalar-tensor action is given by [3] 

1 



SsT — \ dtdrr 



where 



and 



= \ [n<5^. - + n^„n", - nn^ + ^(n^ - n^^n^jj^. 

Also in this case, we consider perturbations around our self-accelerating configuration. 



so that the scalar-tensor action becomes 
where 



o _ o(l) , c.(2) 
i>ST — JsT + "^STi 



S 



1) 
ST 



dtdrr 



fit, r) - ^Sx"" [£ fit, r) r^l^ + 5x^,X<^T 



with 



S 



(2) 
ST 



qS 1 Qint 1 qT 



oS 



s 



ST 



-^Sx^''i£5x)^.., 



and X*^"-* the n"^ perturbation of (5.14). Given our field configurations, the non-vanishing components of the 
tensors X'*' read 

5.22) 
5.23) 





= --5^' 
8 






/J 


r 



2(57r'' 



— Sir 5ti" + 



25tt' 



5.12) 

5.13) 
5.14) 

5.15) 
5.16) 
5.17) 

5.18) 
5.19) 
5.20) 

5.21) 



5.24) 
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Starting from the previous expressions, we can show that 



0(1) , c.(i) _ 



0, 



(5.25) 



as expected, since first order perturbations should vanish using the equations of motion for the background. 
Furthermore, the pure scalar field perturbation simplifies to 



qS _ 



■ dtdr 



(5.26) 



Finally, we only need to focus on the scalar part of the tensor perturbations with respect to the 2 sphere, since 
it is the one that has non-trivial dynamics due to the coupling with the scalar degrees of freedom. It is convenient 
to use the gauge [27] 

/ -Ht -Hi 0\ 
Hi Hr 



The tensor-scalar coupling is then given by 

dtdr 

while the Einstein- Hilbert action reduces to 





V 




/ 



(5.27) 



Qint 



^[ -3 



9 + 16Qg)(57r'i7, 



S'Jy — I drdtr^ 



1 

2^ 



1 

2r 



2rHrHi + ]^Hl - ^Ht{2rHry 



(5.28) 



(5.29) 



By taking variations with respect to Hi and Ht, we find that Hr is not dynamical. Thus, tensor degrees of freedom 
do not play any role in the dynamics of the system in the spherically symmetric case, and their dynamics can be 
consistently neglected. 

5.3 The kinetic terms, and the inevitable ghost degree of freedom 

Collecting the previous results we obtain the second order action for the scalar and vector field perturbations, given 

by 



r dtdr 



6O; 



8 (3+^9 + 16g2)2 



St:' 



(3+ V9 + 16g§)2 



5tt' 



n 



4go 



3 + v/9Tl6Q2 

It is straightforward to diagonalize the kinetic terms, and identify ghost-like degrees of freedom. Defining 

^ ^ 3 ^ 6Qg 2Qo 

8 (3+ v/9-H l6g2)2' P 3 + ^9 + 16QI ' 

we get 

a P \ d"^ f Stt 



S 



(2) 



dtdrr 



{67:, SAt) 



P dt^ \ SAt 



Then by using the the following definitions 

Stt = Sir + SAt, 

SAt 



-OTT 



2/3 



2/3 



SAt 



(5.30) 

(5.31) 
(5.32) 

(5.33) 
(5.34) 
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the kinetic term reads 



5(2) = / dtdrr 



[5^, 5 At) 



Ai 





A2 



5At 



where 



Ai = ^(a + ^^2+4/32), A2 = ^(a- Va' + 4/32). 



(5.35) 



(5.36) 



A ghost is present if one of the eigenvalues Ai is positive. It is easy to see that Ai > 3/8 and A2 < 0, with both 
equahties satisfied when Qq = 0. Thus Stt is always a ghost while 6At is a normal mode, regardless of the size 
of the vector charge Qo- Notice, however, that when Qo = the kinetic term for the mode 6 At vanishes and we 
enter in a regime of strong coupling. The background vector field with Qo ^ then alleviates the strong coupling 
behaviour providing a healthy kinetic term for the vector fluctuations. In the next section, we investigate whether 
these conclusions change when as and 04, and thus higher powers of /C, are included in the potential. 

6 The case of non-vanishing and 

In this section we would like to discuss what happens in the general case, when higher powers of K. are included 
in the massive gravity potential (2.2). The inclusion of such terms changes only the coefficients in front of each 
contributions, but not the general structure of the solutions and the Lagrangian in the decoupling limit. However, 
the existence of a ghost does depend in a subtle and surprising way on these coefficients, as we will see in what 
follows. 

6.1 Lagrangian in the decoupling limit 

In order to construct such Lagrangian using the material developed in Section 4, we first write the potential (2.2) 
in terms of M, resulting in 



U 



+a4 3 



- e) (Vm) + as [24 - 2(Af3/2) + 3(M) ((Tm) - 2) - (Vm) (i8 + ((VM) - e) (%/M)) 
(8 + ((M) - 4) (M) - 8(\/M)) - 6(Af2) + 8{M^^^) ((Vm) - l) + {^) ( - 6(Af) ((\/a7) - 2) 



+ (^/M> (12 + (-4 + (^/M)) (^/M >) ) 



+ 12 - (M) 



(6.1) 



Notice, that not only the traces of \/M and M are needed, as in the a^ ^ a4 = case, but also the trace of 
][,f3/2^ Therefore, it is necessary to solve completely for E using equation (4.22), and not only for its trace, as we 
did in section 4. In the case of the most general spherically symmetric Ansatz (5.1), one gets a similar interaction 
Lagrangian to that of (4.37), but including terms with 03 and a^. Since its expression is long, we refer the reader 
to Appendix C for its explicit form. However, we can mention that the vector-scalar coupling is very similar and 
has the same expression in the denominator as before, i.e. (2 + tt — tt"). Therefore, the same concern about higher 
derivatives applies; on the other hand, the generalised self-accelerating solution is such that higher derivative terms 
vanish when perturbing around it. In what follows we construct such self- accelerating solution. 



6.2 Self-accelerating solution 

The most general static spherically symmetric solution can be constructed in the same way as for the case studied 
in section 3. Again, for the full Lagrangian (2.1), there are two branches of solutions: one with a diagonal metric 
and the other with an off-diagonal component. In contrast to the = a^ = case, the diagonal branch presents 
perturbative solutions which do not decay at infinity, but we will not focus on them, since they do not recover 
general relativity, via the Vainshtein mechanism, for small radius [5, 6, 13, 26]. The non-diagonal branch is very 
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similar to previous case, and a full solution can be constructed with the same structure. Details are given in 
Appendix A. For the purposes of this paper, we only need the decoupling limit of this class of solutions, and which 
is given by (see Appendix B) 



= -A3 

1 

^2 



(1 + 3a3 + 205) + (1 + Sag + a^fA' 



6(l + 3a3 + «5)2 



1 - 



(2 + 3Q3 +Q5) 

(1 + 803 + 05) 



r(A, Qo, as, 04) 



1 



2(l + 3a3+a5) 



A Qo . 2 



where 



r(A, Qo,a3,a4) = 

We have assembled the parameters and a4 into a combination called a^, with 



3(l + 3a3 + a5)^Q§ 

(2 + 3^3 + a5)2[(l + 3^3 + 2a5) + (1 



303 + 05)^ A]' 



1 + 3^3 + 9al - 12a4 



(6.2) 



(6.3) 



(6.4) 



to simplify considerably the expressions. Notice, that this solution includes the possibility of a bare cosmological 
constant A in our formulae, that we added for completeness, and using the normalisation (2.9). Moreover, as 
should be real for the solution to exist, resulting in 



1 + 3a3 + 9al - 12^4 > 0, 



(6.5) 



and it can have both signs, depending on the branch of the square root one considers. Figure 1 shows the parameter 
space in which the solution (6.2) exists for both, positive and negative signs, of as. Therefore, the limit — = 
can only be taken in the positive square root branch, otherwise one is lead to singular expressions. The sign of the 
combination [l + 3a3 + 2as + (1 + 3a3 + a5)^A] sets whether the solution is de Sitter or anti-de Sitter, and one 
can associate an effective cosmological constant^ given by 



A, 



eff,± 



A 



1 + 3a3 ± 2^1 + 3a3 



■9al 



YlaA 



(6.6) 



(1 + 3a3 ± v^l + 3a3 + 9a\ - Yla^Y ' 
In the AdS case, there is an extra constraint on Qo due to the argument in the square root of (6.3), and given by 

(2 + 3a3 + as)^ (1 + 3a3 + 2as) + (1 + 3a3 + a^f A 



< 



3(l + 3a3 + a5)4 



(6.7) 



Although we derived these decoupling solutions from the exact spherically symmetric solution, other exact 
solutions in the full theory may have the same decoupling limit. For example, the open FRW solution found in 
Ref. [10] also reduces to Eq. (6.2) with Qo = in the decoupling limit. 

As expected, the previous expressions (6.2) solve the equations of motion obtained from the decoupling limit 
Lagrangian, which was discussed in the previous Section and whose full expression is given by (C.2). Now, we would 
like to understand the stability of these, more general, self-accelerating solutions, which include vector degrees of 
freedom and a bare cosmological constant. 



6.3 Linear perturbations: choosing between a ghost or strong coupling 

We now study spherically symmetric perturbations around the background (6.2) as we did for the = 04 = 
case in Section 5. The procedure is exactly the same, and we refer the reader to Appendix C where the explicit 

®In the expression for the effective cosmological constant we have not dropped the normalisation of (2.9). 
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Positive Case 



Negative Case 



a. 



Ill dS 
E AdS 



1 1 - Brane 



04 



^ III dS 

E AdS 



a,. 



Brane 



»,1 04 



Figure 1: The dS and AdS regions of the self-accelerating sohition in the decouphng hmit (6.2). Left (right) plot 
is for 015 > (a5 < 0). The white exclusion area is associated to the positiveness of the argument in the square 
root of Qfs. In each plot, half of the line + 3a4 = has to be removed to avoid the curvature of the space-time 
from diverging, with the middle point being that of the brane action found in Ref. [6]. 



calculations are shown. Here, we only present the results. The Lagrangian controlling the perturbations is 

(2) _ I 2,,, / [1 + 3q3 + 2^5 + (1 + 3^3 + «5)^A] [Stt'^ ~ r(A, Qo, a3, a4)STr'^] 
^ - (l + 3a3 + a5)2[l + r(A,Qo,«3,a4)] ^^'^^ 



(2 + 3a3 + a5)[l + r(A, Qo, "3, "4)] 

where T is defined in (6.3). Therefore, we only need to consider fluctuations in Stt, and in SAt for the vector. Their 
kinetic terms can be straightforwardly diagonalised, and the resulting action is again (5.35), but with eigenvalues 

X12 = |(q^ + ± 4/3^), where a and /? are now given by 

as Tl -I- 3a3 + 2a^ + (1 + 3a3 + a5)^Al r(A, Qq, as, 0:4) 
a = 



(1 + 3a3 + a5)2[l + r(A, Qo,a3,a4)] 



,1/2 



as ([2 + 3a3 + a^]^ [r(A, Qp, ag, a4)^ - l] [l + 3«3 + + (1 + 3a3 + a5)^A] ) 

v^(l + 3a3 + as)[l + r(A,Qo , a3,a4jj 

In order to investigate the existence of ghost degrees of freedom, we have to study the sign of A^; if at least one of 
the two eigenvalues is positive, a ghost is present. In order to study the ghost issue, it is convenient to focus on 
the product of the eigenvalues, whose expression reduces to 

aj [1 + 3a3 -I- 2as + (1 + 3a3 + as)^A] [r(A, Qp, as, aj) - 1] 

3(l + 3a3 + as)2[l + r(A,Qo,a3,a4)] ^' ' 

It is simple to convince oneself that the previous product is always negative when Qo 7^ 0; thus a ghost is always 
present around the self-accelerating solution with the non-vanishing vector field because one of the modes have 
kinetic terms with the wrong sign. The result is independent on the values of a^ and A, as long as we have dS (or 
AdS) symmetry. On the other hand, the kinetic term for the other mode is positive definite and always healthy. 

On the other hand, when Qq = (F = 1), the product (6.10) vanishes. The ghost can be avoided but a strong 
coupling appears since the kinetic term for one of the fluctuations goes to zero in this limit. It is not the purpose 
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of this paper to investigate in detail the physical consequences of this strong couphng regime. Here, we are only 
interested in the conditions that avoid the ghost, generahzing the analysis of [8] to the case of non- vanishing bare 
cosmological constant. Therefore, when Qo = the eigenvalues read 

Ai,2 = 4^(^i + + ar)^ [05(1 + Sag + 2^5 + (1 + Sag + a5)2A) ± jasHl + 803 + 2^5 + (1 + Sag + a5)^A|] . 

(6.11) 

We assume the Hubble parameter squared is positive (de Sitter regime), which from (6.2) translates into the 
condition 

l + 3a3 + 2a5 + (l + 3a3+a5)^A>0. (6.12) 

Comparing this with Eq. (6.11), we conclude that the only way to avoid the ghost is to choose as < 0. To 
summarize, the conditions to avoid the ghost in a de Sitter regime are 

— 1 — Sag + 2|at;| 

A>— — and a5<0. 6.13 

(1 + 3a3 - \a5\y 

These conditions generalize the results previously obtained in the literature, and when A = they reduce to those 
in [8]. 



7 Conclusions 

In this paper, we presented new self-accelerating solutions for a family of non-linear massive gravity models in the 
decoupling limit. The dynamics of vector degrees of freedom of the massive graviton is fully taken into account, 
playing an important role in characterizing these configurations. We found two branches of solutions, which describe 
a de-Sitter spacetime in different regions of the parameter space (a3, a^). In order to study the dynamics of linear 
fluctuations, we developed the tools to obtain a general Lagrangian in the decoupling limit, which includes vector 
degrees of freedom. The resulting Lagrangian contains an infinite number of terms, which lead to non-polynomial 
interactions between the vector and scalar degrees of freedom. 

We then studied linear fluctuations around these self-accelerating configurations, finding two different be- 
haviours. When the background vector field is switched on, a ghost mode is always present in the spectrum 
of propagating modes, regardless of the choice of parameters characterizing the theory. In contrast, when the 
background vector field is set to zero, our configurations reduce to those analysed in [8]. In this case, ghosts may 
be absent in one of the two branches of solutions and in some regions of parameter space. This ghost-free branch 
does not exist for the simplest case with a3 = a4 = 0. However, the coefficient in front of the kinetic terms for the 
vector fluctuations vanishes, leading to a strong coupling effect. 

A bare cosmological constant can be fully included in these self-accelerating configurations, and the same 
conclusions hold. There is ghost mode if a background vector field exists, and the ghost may be removed if this 
background field is not present. The avoidance of the ghost is again only in one of the two branches of solutions, 
and for a positive bare cosmological constant the parameter space with a ghost-free model is enlarged. 

Our analysis can be extended in many interesting directions. It is important to investigate in more details 
the Lagrangian in the decoupling limit with vectors. It would be interesting to analyse how matter couples to 
the degrees of freedom controlled by this Lagrangian and understand the role of vectors. Regarding the self- 
accelerating solutions, it would also be interesting to understand whether symmetries exist that automatically set 
the background vector field to zero and avoid ghost modes. Imposing Lorentz invariance on the solution of the 
scalar mode, tt, can achieve this goal; but a more comprehensive analysis of possible symmetries is needed. It 
is also crucial to clarify the physical consequence of the strong coupling regime in the vector fluctuations, that 
one encounters when the vector background profile is switched off. As stated in [8], quantum corrections might 
be able to induce a healthy kinetic term for the vector; it would be interesting to understand this effect in more 
detail, and to discuss its consequences for the self-accelerating configurations. Finally, it is encouraging that there 
exist ghost- free self-accelerating solutions in the decoupling limit, but this does not necessarily ensure that the 
full non-linear solution is stable, especially when a bare cosmological constant is included. This last point may 
be relevant for inflationary cosmology within this theory of massive gravity. A detailed analysis of full non-linear 
self-accelerating solutions away from the decoupling limit will be presented in a future publication [28]. 
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A General exact solution 



From the general Lagrangian (2.1), and using the non-diagonal ansatz (3.1) together with Einstein equations (3.3), 
one can show that there are two branches of solutions for non- vanishing and a^. These solutions were studied 
in detail in [6], however, here we only consider the branch with a non-diagonal metric, where analytic solutions 
can be found. Since the combination -y/l -I- Saa -|- 9a| — 12a4 is always present in the solution of this branch, it 
is convenient to map the (03,0:4) parameters into (03,05), where 05 = 1 -I- 3o3 + 9o| — I204. In this new set of 
parameters, the combination D{r) Gtt + C{r) Gtr — 0, fixes i? as a function of r in the following way 



B{r) = bor' = 



2 _ (1 + 303 05)2 2 



(2 + 303 + 05)2 



(A.l) 



The rest of Einstein equations give 



^W = ^(l-P), A(r) = ^(p + 7 + l), D{r) = ^Ao-A{r)C{r), 



(A.2) 



where 



P = 



3(2- 



3o3 + 205) 
- 3o3 + O5)' 



2 2 

-m r , 



7 + 1 = 



(1 + 303 + 05)4 
Ao(2 + 303+05)4 



(A.3) 



Just like in the 03 = 04 (05 — 1) case, there are two integration constants, c and Aq, but in order to have a 



positive argument for the square root in D{r), Aq has to run from Aq 



to AJ)"°^ 



hi. In this paper, we focus 



on the massless case c = only, which describes the static patch of the de Sitter or Anti-de Sitter spacetime. 



B The decoupling limit 

If one naively takes the decoupling limit (2.11) of the static configuration (3.6) [or more generally of (A.1)-(A.2)], 
one gets a divergent result. This is due to scalar and vector contributions which would be zero in this nai've limit. 
Therefore, in order to extract a finite answer, one should look for a coordinate transformation that, after expanding 
h^v = — ill terms of m, leaves the first non- vanishing coefficient to be that of m?, so that the limit to — 
is well defined. For our case, this transformation is precisely given by 

{t,r)^ f^AAi + f , , (B.l) 



\ "'0 

where h^, D and C are given in (3.6) [or more generally in (A.1)-(A.2)]. The metric field h^n, after being canonically 
normalised using the definition (2.9), reads 

_ (1 + 303 + 205) 2 
^''"3(1 + 303 + 05)2'^^^" ^ ' 

Moreover, the transformation (B.l) induces the following components of the Stiickelberg field tt'', 
_ ( I + 303 + 05 1 /(I + 303 + 205) [(1 + 303 + 05)4 - (2 + 303 + 05)4Ao] . 2, 



V(2 + 3o3 + 05)VA^ J 2\ 3Ao(l + 303+05)4(2 + 303 + 05)2 ^ ' 

— -1— r + 0(TO2), (B.3) 
1 + 303 + 05 
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which in turn, can be decompose into scalar and vector modes. After being canonically normaUsed using (2.9), the 
vector and scalar modes reduce to 



2(l + 3a3 + a5) 



(2 + 3a3 + a5)Jl 



3(l + 3a3 + a5)^Qg 
(2 + 3a3+a5)2(l + 3a3 + 2a5)' 



Qo 



(B.4) 



where we have decided to measure the vector charge using Qo — a real positive number — instead of Ag (which 
was bounded from to 5q). The following equation relates Qo Aq, 



Ao = 



(1 + 3a3 + q;5)^(1 + 3q;3 + 2Q5) 

3Q2(2 + 3^3 + "5)^(1 + 3^3 + as)" + (2 + 3^3 + a:,)^{l + 3^3 + Sag) ' 



(B.5) 



Finally, one can take the decoupling limit of expressions (B.2) and (B.4), and get a well defined answer. 

One can induce further coordinate transformations, of order m^, which take the metric to any desired form in 
the decoupling limit and consequently do not affect the vector nor the scalar. Using this freedom, one can write 
the canonically normalised metric field ft.^^ in a covariant conformal form, namely 



(l + 3a3 + 2a5)A3(r^-t^) 
6(l + 3a3+a5)2 



+ 0{m)] rj^^. 



(B.6) 



Therefore for (1 + 3a3 + 205) positive (negative) one gets de Sitter (Anti-de Sitter). In the case of AdS, one has 
an extra constraint on Qq, from the square root argument in (B.4), given by 



< 



(2 + 3q;3 + a5)2|l + 3^3 + 2a5 
3(1 + 3^3 + 05)4 



(B.7) 



One can include a bare cosmological constant, and still get a solution for the non-diagonal branch with a very 
similar form as (A.1)-(A.2). The details of this solution were described in [6], but here we will only describe its 
decoupling limit. Therefore, the non-diagonal branch solution for the action (2.1), which includes a cosmological 
constant A (normalised as in (2.9)), has the following decoupling limit 



-A3 
1 / 



(1 + 3a3 + 2^5) + (1 + 3a3 + a5)'^A 
6(1 + 3(13 + 05)2 



(r2 - t^)7^^,, 



1 



(2 + 3^3 + as 
(1 + 3a3 + as) 
1 




3(l + 3a3 + as)'iQ2 



(2 + 3a3 + as)2[(l + 303 + 2az) + (1 + 3a3 + a5)2A] 



2(l + 3a3 + a5 



-A3r^ 



Ao = 



Qo 



In this case, the relationship (B.5) between Aq and Qq, gets modify to 

(1 + 3a3 + a5)4[(l + 3a3 + 2as) + (1 + 3a3 + asfA] 



Ao = 



3Q2(2 + 3a3 + as)2(l + 3a3 + as)* + (2 + 3a3 + as)4[(l + 3a3 + 205) + (1 + 3a3 + as)2A] 
and the AdS bound (B.7) to 

(2 + 3a3 + as)2 (1 + 3a3 + 2as) + (1 + 3a3 + a5)^A 



(B., 



(B.9) 



3(1 + 3a3 + 05)1 



(B.IO) 
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C Decoupling limit Lagrangian with spherical symmetry 

Let us consider the following Ansatz 

h^..^ f{t,r)7j^„ TT^TT{t,r), = {At{t,r), Ar{t,r), 0,0). (C.l) 

By inserting this into the action (2.1), one can take the decouphng hmit to obtain the following Lagrangian 



^dec 



+ C'^'if) + F2in)fit, r) + F3(7r, A, Ar) 



(2) 

where we have set A3 to one. The first term, Fi, is a total derivative and Cp^ {f) is the quadratic expansion of 
Einstein piece with a bare cosniological constant A, which reduces to 



(C.2) 



4'H/)-|:(/"-/')-4A/. 



(C.3) 



The term with F2 is simply given by [3] 



F2 - ^Tr(AW) = -3('if-7r"-2-) + [(l + 3a3f-ai] 



i=l 



2(7r'2 - ifTr") -{^- tt" 



r 



+4(1 + 3a3) 



— ( 27r" + — I - TT I TT 

r \ r 



(C.4) 



Finally, the interaction term between the vector and the scalar, can be obtained using the techniques described in 
section 4, but solving for E, and not only tiE, since M'^/^ is also needed when cubic or higher powers of K. are 
included. It is an straightforward generalisation to the = = case, and the final results is 



2 + TT - tt" 



1 + ^ - tt" + -(-1 + 303(1 + n- n")y 



[3^3 - (1 + 3^3 + Qal -al){l + n- tt") 



^'21 



Ar ~1~ A^ 



n 
r 



1 + — ( 2 + 6a3 + ^ (1 + 6a3 + 9a^ - al] 



-2A'rAt {2 + n- tt" 



1 + ^ ( 6^3 + ^ (1 + 3^3 + 9al - al) 



-Ar\ a; -1 + 3^3 if I + (1 + 3a3 + - a^) 



-{Ar + A;)7t' 3^3 + — (1 + 3^3 + 9^3 - ai^ 



+it 1 + (1 + 3a3 + Qaj - al 



3q!3 tt 



2 



1 + 3a3 {-n + tt") + (l + 3^3 + 9al - al) (tt'^ - nn") 



At 



(C.5) 



One can check that the self- accelerating configuration (6.2) solves the equations of motion for this Lagrangian. 



C.l Perturbations 

In order to perform the perturbation analysis, we define the following perturbations 

h).^ ^ f{t,r)or]^^ + 6xt.^, tt = ttq + <57r(t, r), A^ ^ {Ato{r) + SAt{t,r),dAr{t,r), 0,0), (C.6) 
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around the self-accelerating solution {fo,T^o,^to) given by (6.2). One can decompose the tensor perturbations £_^,y 
in the same way as before (see Eq. (5.27)), and show that do not contribute to the dynamics. The pure vector part 
is a total derivative as before, and only the vector-scalar mixing and pure scalar part play a role. The final action 
to second order in perturbations is 



^(1) 



S 



vs 



(C.7) 



where subscript represent the origin of such term (ST=scalar-tensor, SV=scalar-vector) and the upper label refers 
to pure scalar (S) or vector-scalar (int). The linear perturbations 



'ts 



as [(1 -t- 3a3 



(2 + 3^3 + a5)2yA^] [1 + Sag + 2a5 + (1 + Bag + a5)^A] Sn'r^ 



3(1 + 3a3 + a5)3(2 + Sa^ + a^)^^ 



(C.8) 



as [(1 + 3a3 -t- as)^ - (2 + 3a3 + as)^^^;;] {rSn" + 26tt'){1 + Sag + 2a5 + (1 + 3a3 + a5)'^A)r{r^ - t^) 



6(1 + 3a3 + as)3(2 + 3a3 -I- a5)v^ 
cancel after an integration by parts, as expected. The quadratic scalar parts reads 

as (l + 3a3 + 2as + (Sag + a^ + ifk) 



Syg + S^g 



r dtdr- 



/A(^(3a3 4- as -I- 2)2(3a3 + + 1)2 + (3a3 + + 1)^ 



/Ao(3a3 -I- as -I- 2f5TT'' 



-(3a3 + as + If Sir' 



(C,9) 



meanwhile the vector-scalar mixing term is 

2a5 ^((1 + 3a3 + as)^ - (2 + 3a3 + as)4Ao)(l + Sag + 2as + {I + 3a3 + as)2A)r2 



Qmix 

Jvs — 



- 5 At Sn 



V^il + 3a3 + as)((l + 3a3 + as)^ + (2 + 3a3 + a^fV^o 

(CIO) 

Therefore, the action in the main text (6.8) is obtained from the last two contributions, but expressed in terms of 
Qo instead of Aq, which is achieved by equation (B.9). 
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